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Nomenclature,  States  and  Inputs 


State 

Input 

vt 

Vehicle  velocity 

<5e 

Control  surface  deflection 

a 

Angle  of  attack 

<f> 

Fuel  to  air  ratio 

Q 

Pitch  rate 

Ad 

Diffuser  area  ratio 

h 

Altitude 

Cowl  lip  position 

e 

Pitch  angle 

Vi 

ith  Generalized  elastic  coordinate 

Vi 

Time  derivative  of  r]i 

Table  1.  States  and  inputs  for  the  model 


I.  Introduction 

Air-breathing  hypersonic  vehicles  may  eventually  allow  dramatic  reductions  in  flight  times  for  both 
commercial  and  military  applications.  Direct  access  to  earth  orbit  without  the  use  of  separate  boosting  stages 
may  also  become  possible  as  scramjet  powered  aircraft  enter  service.  While  numerous  challenges  remain, 
recent  successes  with  NASA’s  X-43a  and  renewed  research  activities  throughout  the  aerospace  community 
suggest  that  this  technology  is  well  on  its  way  to  assuming  a  role  in  the  next  generation  of  aviation. 

The  design  of  controllers  for  air-breathing  hypersonic  vehicles  requires  the  controls  engineer  to  overcome 
strong  couplings  between  the  propulsive  and  aerodynamic  effects  while  also  addressing  the  significant  flex¬ 
ibility  associated  with  the  slender  geometries  required  for  these  aircraft.2, 14,15  As  a  result,  all  the  control 
schemes  proposed  to  date  have  been  limited  in  their  scope  to  the  longitudinal  dynamics  of  hypersonic  vehi¬ 
cles.  A  wide  range  of  control  laws  have  been  developed  for  linearized  versions  of  hypersonic  models,1,12,14-16 
while  a  few  have  also  attempted  to  incorporate  guidance.5, 17  Tournes  et  al.  employed  a  nonlinear  variable 
structure  control  approach  in  Ref.  18,  while  several  other  nonlinear  control  approaches  have  been  completed 
on  nearly  non-flexible  models.6,11,19,20 

This  paper  deals  with  nonlinear  control  of  the  model  of  a  flexible  air-breathing  hypersonic  vehicle  de¬ 
veloped  by  Bolender  and  Doman,2  which  employs  compressible  flow  theory  to  obtain  a  more  complex  and 
potentially  more  complete  model  than  those  developed  in  earlier,  similar  efforts.3,4  Following  the  successful 
design  of  linear  controllers8  for  the  model  considered  in  this  paper,  the  next  logical  step  is  the  synthesis  of  a 
nonlinear  control  law  for  this  novel  plant.  In  this  initial  study,  a  traditional  feedback  linearization  approach 
will  be  pursued. 

Intricate  interactions  between  the  structural,  aerodynamic,  and  propulsion  system  equations  result  in 
exceedingly  complex  expressions  for  the  aerodynamic  forces  and  coupling  terms  between  the  flexible  states 
and  the  other  equations  of  motion.  When  written  in  a  form  suitable  for  control,  these  equations  are  clearly 
analytically  untractable,  even  with  the  aid  of  a  computer  algebra  system.  The  intractable  model  will  be 
referred  to  as  the  truth  model  and  will  be  used  solely  for  simulation  and  control  design  verification. 

The  first  step  in  the  control  design  process  is  the  creation  of  a  nominal  model  which  approximates  the 
behavior  of  the  truth  model  with  reduced  complexity.  The  outputs  to  be  controlled  are  chosen  as  the  vehicle 
velocity  and  altitude.  The  control  inputs  are  the  elevator  deflection  and  fuel  to  air  ratio.  The  resulting 
non-linear  nominal  model  is  at  first  glance  very  similar  to  those  given  in  previous  works.19,20  However,  as 
opposed  to  the  models  considered  in  Refs.  19,20,  the  addition  of  flexible  states,  increased  complexity  of  the 
engine  model,  and  additional  parasitic  couplings  between  the  control  effectors  and  the  aerodynamic  forces 
create  an  unstable  zero  dynamics  and  drastically  increase  the  complexity  of  the  vehicle’s  governing  equations. 

Rather  than  attempting  to  design  control  directly  for  this  nonminimum-phase  system,  the  approach 
of  approximate  feedback  linearization  is  employed.9  The  flexible  states  and  a  set  of  weak  couplings  are 
removed  from  the  nominal  model  to  create  a  simplified  nominal  model  which  can  be  shown  to  posses  full 
vector  relative  degree  with  appropriate  dynamic  extension.  Thus,  a  straight-forward  design  of  a  feedback 
linearization  control  law  can  be  completed  for  this  simplified  nominal  model.  An  outer  loop  controller  is  then 
designed  using  linear  quadratic  regulator  (LQR)  methods  with  integral  augmentation  and  model  reference 
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control.  Simulation  studies  reveal  that,  despite  the  simplified  model  of  the  plant  dynamics,  this  control  law 
performs  admirably  on  both  the  nominal  and  truth  models. 

The  paper  is  organized  as  follows:  section  II  describes  the  construction  of  the  nominal  model  and  the 
system  identification  process  that  is  employed  to  match  it  to  truth  model  data.  In  section  III,  the  simplified 
nominal  model  is  presented.  The  control  law  design  is  described  in  section  IV,  and  simulation  results  for 
the  truth  model  are  given  in  section  V  along  with  simulation  results  for  small  plant  parameter  variations. 
Finally,  an  alternative  choice  of  outputs  and  the  corresponding  simulation  results  are  given  in  section  VI, 
while  section  VII  concludes  the  document  with  a  brief  summary  of  the  results. 


II.  Nominal  Model 

In  this  section,  we  derive  a  control-oriented  version  of  the  model  of  the  longitudinal  dynamics  of  an  air- 
breathing  hypersonic  vehicle  developed  in  Ref.  2.  The  simplification  is  necessary  to  overcome  the  analytical 
intractability  of  the  original  model,  which  will  be  henceforth  referred  to  as  the  truth  model.  For  the  sake  of 
simplicity,  the  control  oriented  model  will  be  referred  to  as  the  nominal  model.  After  the  synthesis  of  the 
model  is  completed,  its  suitability  for  control  design  is  evaluated. 


A.  Model  Development 

The  nonlinear  longitudinal  equations  of  motion  for  the  model  developed  by  Bolender  and  Doman2  are  given 
as 


Vt 


a 


lyyQ 


h 

0 


kiVi 


—  (T  cos  a  —  D)  —  g  sm(6  —  a) 
m 

— 77  (-T  sin  a  -  L)  +  Q  +  cos (0  -  a) 

mVt  Vt 

M  +  -01771  + 

Vt  sin (0  —  a) 

Q 


2  ,  AT  7  M 

-2Cicji?7i  -  0^7/1  +  Aq  -  - - - - 

J-yy  J-yy 

0/.  2  ,  AT  7  M  Mi  m 

-2C,2^2m  ~  u2ri2  +  N2  -  W2j - j — , 

yy  * yy 


(i) 


(2) 


where  the  system  state  xn  G  M9  is  comprised  of  the  nine  state  variables  Vt,  o,  Q,  h,  61,  771 ,  772,  and 
1)2  which  are  defined  in  Table  1  along  with  the  four  control  inputs.  This  formulation  of  the  system  differs 
slightly  from  that  given  in  Ref.  2  in  that  the  temperature  change  across  the  combustor  input  ATo  has  been 
replaced  with  the  functionally  equivalent  control  input  <f>,  the  fuel  to  air  ratio.  The  change  was  made  since 
<F  is  a  more  natural  input  choice  for  aerospace  design.  The  outputs  to  be  controlled  are  Vt  and  h.  Since  only 
two  control  inputs  are  required  for  the  feedback  linearization  design,  Xd  and  Ad  will  be  fixed  at  the  nominal 
values  Xd  =  0  and  Ad  =  1  for  the  remainder  of  this  study.  Table  2  defines  several  physical  constants  derived 
from  the  vehicle  geometry,  aerodynamic  conditions,  and  assumed  elastic  mode  shapes.2 


Physical  Constant 

m 

Vehicle  mass 

lyy 

Moment  of  intertia 

Constrained  beam  coupling  constant  for  r]i 

h 

1  + f- 

1yy 

Ci 

Damping  ratio  for  elastic  mode  rji 

Uli 

Natural  frequency  for  elastic  mode  rji 

Table  2.  Physical  constants  included  in  the  equations  of  motion  (1)  to  (3) 
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The  lift  L,  drag  D ,  thrust  T,  pitching  moment  M,  and  the  two  generalized  forces  TVi  and  A2  are  complex 
algebraic  functions  of  both  the  system  state  and  the  inputs  that  must  be  simplified  to  render  the  model 
analytically  tractable.  While  attempts  to  represent  these  force  mappings  with  cubic-spline  interpolated 
lookup  tables  were  very  successful,7  these  approximations  are  not  suitable  for  analytical  analysis  with  the 
tools  of  feedback  linearization.  Thus,  drawing  on  earlier  work  done  with  a  non- flexible  vehicle,19,20  the 
following  approximations  to  these  mappings  will  be  adopted  for  the  nominal  model 

L  «  \pV2SCL  (3) 

D  «  \PV?SCD 

M  PS  ZtT  +  -pV2  Sc[CM,a  +  Cm,c$J 

Ni  ps  Nf  a2  +  N?a  +  N° 

N2  ps  Nfa2+N°a  +  N^5e  +  N$. 


The  expressions  for  L  and  D  are  the  same  as  those  given  in  Ref.  20.  The  expression  for  M  contains  the 
additional  term  ZtT  to  account  for  the  pitching  moment  produced  by  the  underslung  scramjet  engine  in  the 
model.2  Finally,  the  numerical  studies  completed  by  Groves  and  Sigthorsson7  support  approximating  the 
Ni  force  as  a  quadratic  of  the  angle  of  attack  and  the  N2  force  as  a  quadratic  of  the  angle  of  attack  with 
a  linear  elevator  deflection  term.  The  Q  coupling  term  to  the  moment  has  also  been  neglected,  because  it 
proved  inconsequential  for  this  vehicle.  The  C*  coefficients  are  given  as  follows 


p  = 

~(h  -  h0) 

Po  exp(  ^  ) 

Cl  = 

Cta  +  C^Se  +  Cl 

Cd  = 

Ch  0 ?  +  C^jQ!  +  C£  S2  +  CpSe  +  Cj) 

Cm, a  — 

^M,aa2  +  C M,aa  +  ^M,a 

Cm, be  = 

ce(6e  -  a). 

The  expressions  for  the  C*  coefficients  are  similar  to  those  in  Ref.  20,  with  minor  notation  changes  and 
the  addition  of  elevator  terms  in  the  lift  and  drag  coefficients.  A  simple  exponential  model  for  the  air  density 

p  is  also  included. 

The  scramjet  engine  included  in  the  model  produces  a  thrust  force  T  which  depends  strongly  on  the 
states  h,  Id,  and  a,  along  with  the  input  <f>.  The  mapping  is  approximately  cubic  in  the  angle  of  attack, 
while  each  coefficient  of  this  polynomial  is  a  linear  function  of  <F.  The  following  approximation  for  T  is 

adopted  for  the  nominal  model 

T 

«  Cfa3  +  Cfa2  +  C%a  +  C £ 

rioP 

Lyrji 

=  fa  4>  +  fa 

/^ia2 

=  fa&  +  fa 

/^iOi 

Lsrj1 

=  fa$  +  fa 

=  +  fa 

Q 

=  \pV?.  (4) 

The  eight  fa  coefficients  vary  with  the  dynamic  pressure  q  and  h.  Since  these  values  change  on  a  much 
slower  timescale  than  4>  and  <a,  the  values  of  fa  are  assumed  to  be  constant  for  control  design.  In  imple¬ 
mentation,  the  coefficients  are  obtained  using  a  cubic-spline  interpolated  lookup  table  based  on  coefficients 
derived  for  400  different  flight  conditions. 

All  of  the  force  approximations  given  in  equations  (3)  to  (4)  are  linear  in  their  fit  parameters.  The  values 
for  po,  hs ,  ho,  zt ,  S ,  and  c  are  known  constants.  The  other  parameters  have  been  obtained  through  a  least 
squares  approach  by  fitting  the  functions  to  truth  model  data  calculated  over  the  following  ranges 


h  e  [85000, 135000] 
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$  G  [0.1, 1.2] 

Se  e  [-15,15] 
a  G  [-10,10] 
q  G  [500,2000]. 

Naturally  the  variation  in  q  is  obtained  by  appropriate  variation  of  Vt.  Even  if  a  given  force  approximation 
does  not  depend  on  a  particular  state  or  input,  that  value  will  still  be  varied  so  as  to  capture  its  average 
effect.  The  flexible  states  rfo  and  rji  are  set  to  zero  during  the  calculation  of  the  data  for  the  least  squares 
fits,  since  none  of  the  force  approximations  depend  on  them.  Many  of  the  conditions  on  these  ranges  will 
not  be  physically  realizable  for  the  model.  These  points  are  simply  ignored  by  the  curve  fitting  algorithm. 

The  resulting  curve  fits  use  3.2  million  data  points  and  capture  at  least  the  qualitative  behavior  of  the 
model  over  the  entire  operating  range.  Figure  1  gives  sample  results  for  the  thrust  fit.  The  identification  of 
these  parameters  completes  the  construction  of  the  nominal  model  given  by  equations  (1)  to  (4). 


(a)  Variation  with  respect  to  a  and  h 


(b)  Variation  with  respect  to  a  and 


(c)  Variation  with  respect  to  a  and  q  (d)  Variation  with  respect  to  a  and  Se 


Figure  1.  Thrust  Curve  Fitting  Result  for  a  typical  trim  condition.  Notice  that  6e  has  minimal  effect  on  the 
thrust,  as  expected. 


B.  Relative  Degree  of  the  Nominal  Model 

In  this  subsection,  we  will  briefly  review  the  basic  definitions  associated  with  the  process  of  feedback  lineariza¬ 
tion  and  then  discuss  the  applicability  of  the  technique  to  the  nominal  model.  This  analysis  will  highlight 
the  benefits  of  approximating  the  nominal  model  to  produce  the  simplified  nominal  model  discussed  in  the 
subsequent  section. 
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As  discussed  in  a  variety  of  sources,10,13  consider  a  generic  multi-input  multi-output  (MIMO)  system  in 
which  the  input  and  output  have  the  same  dimension 


m 

x  =  /(*)*£»«« 

i= 1 

Vi  =  hi(x)  (i  =  1  ...to), 

with  x  G  Mn,  u  G  Mm,  and  y  G  Mm.  The  vector  fields  /  and  ^  along  with  the  functions  hi  are  assumed  to  be 
smooth  in  their  arguments,  allowing  arbitrary  derivatives  of  each  to  be  calculated  without  concern  for  their 
existence.  The  following  standard  notation  for  the  Lie  derivative  of  the  function  h  along  the  vector  field  / 
will  be  adopted 


Lfh  =  g/M. 

The  gradient  of  the  scalar  function  h  is  interpreted  as  a  row  vector.  Higher  order  Lie  derivatives,  Lkh  with 
k  G  N,  are  obtained  recursively  using  this  definition.  For  clarity,  the  notation  Lu.h  will  sometimes  be  used 
to  represent  the  Lie  derivative  of  the  function  h  along  the  vector  field  corresponding  to  the  input  U{. 

The  analysis  begins  with  taking  successive  time  derivatives  of  each  output.  Define  r*  such  that  is  the 
first  derivative  of  yi  that  depends  explicitly  on  one  or  more  of  the  inputs  Ui.  This  derivative  can  be  expressed 
for  each  i  as 


V. P  =  Lrfihi  +  Y,La^Lrr1hi)uk. 

k= 1 


The  decoupling  matrix  A(pc)  :  Mn  —>  Mmxm  contains  the  m2  Lie  derivatives  which  multiply  the  inputs  in  the 
m  derivative  expressions  y\l .  Specifically,  A  is  given  as 


A(x)  = 


LgM'rlhi 


L9lUr^K 


L9rnLy-lK 


If  the  matrix  A{pc o)  is  nonsingular  at  xo,  then  the  system  (5)  is  said  to  have  vector  relative  degree  r  at  xo 
where 


r  = 


If  r  =  n  then  the  system  is  said  to  have  full  vector  relative  degree,  and  the  input/output  map  can  be  rendered 
a  set  of  m  decoupled  chains  of  integrators  using  the  feedback  law 


u 


-A 


-l 


Ljhx 


T  T  m 

L  ^ f  '^r 


+  A_1v 


The  new  outer  loop  control  v  G  Mm  contains  the  inputs  to  the  decoupled  integrator  chains.  If  r  <  n  then 
this  process  can  still  be  carried  out,  but  n  —  r  states  of  the  system  will  be  rendered  unobservable  from  the 
output  and  constitute  an  internal  dynamics  for  the  system  under  the  input/output  linearizing  control  law. 
An  extensive  discussion  of  these  results  from  a  differential  geometric  point  of  view  can  be  found  in  Ref.  10. 

Now,  consider  the  particular  case  of  (5)  representing  the  nominal  model  of  the  air-breathing  hypersonic 
vehicle.  Recall  that  the  nominal  model  state  xn  G  M9  contains  the  state  variables  given  in  table  1.  The 
outputs  were  selected  as  y\  =  Vt  and  y^  —  h.  Finally,  the  two  control  inputs  are  Se  and  <F.  For  this 
input/output  combination,  the  matrix  A(xn)  is  singular  over  the  entire  flight  envelope.  In  order  to  obtain 
vector  relative  degree  for  this  case,  the  following  dynamic  extensions13  are  necessary 
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(5) 

(6) 

where  we  have  defined  a  new  input  to  the  system  as 


^  —  — 2^36^3^  —  UJ 3<h  ~b 

5e  =  —2()4Ld4:5e  —  <J\  Se  +  <jj\  Je?c, 


The  values  for  (3,  £4,  CJ3,  and  CJ4  are  selected  to  provide  reasonable  second  order  responses  to  these  extended 
dynamics. 

Considering  this  extended  system  with  state 


%na 


<f> 

i> 

<5e 

Se  . 


the  new  matrix  of  Lie  derivatives  A(xna)  is  nonsingular  over  the  operating  range  of  interest.  Thus,  the 
nominal  model  dynamically  extended  with  (5)  and  (6)  has  a  well  defined  vector  relative  degree.  However, 
since  r  =  7,  the  nominal  model  under  the  feedback  linearizing  control  law  (5)  will  contain  a  six-dimensional 
unobservable  internal  dynamics.  The  following  diffeomorphism  is  used  to  write  this  system  in  normal  form 


Vt 

LfVt 

L)Vt 

h 

Lfh 

L2fh 

L)h 

m 

m 

m 

m 

0 

.  Q  _ 


Symbolic  calculations  in  Mat  lab  verify  that  ^  is  nonsingular  over  the  operating  envelope.  The  four 
flexible  states  rji  and  f]i  along  with  the  pitch  dynamics  6  and  Q  become  an  unstable  internal  dynamics  for 
the  plant.  While  explicit  calculation  of  the  normal  form  in  coordinates  proved  to  be  intractable,  Jacobian 
linearization  about  a  trim  condition  verifies  that  the  internal  dynamics  are  unstable.  It  should  be  emphasized 
that  the  pitch  dynamics  become  unstable  under  feedback  linearization,  even  if  the  flexible  effects  are  removed 
from  the  system.  Since  the  presence  of  unstable  internal  dynamics  substantially  complicates  the  controller 
design,  an  approach  for  removing  the  internal  dynamics  from  the  model  used  for  control  design  will  be 
adopted  in  the  sequel. 


III.  Simplified  Nominal  Model 

The  six  dimensional,  unstable  internal  dynamics  makes  the  nominal  model  unsuitable  for  feedback  lin¬ 
earization  control.  To  avoid  dealing  with  the  internal  dynamics  directly,  a  simplified  nominal  model  which 
enjoys  the  property  of  full  vector  relative  degree  will  be  derived. 
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Despite  their  complex  coupling  to  the  rest  of  the  system,  the  flexible  modes  are  asymptotically  stable 
with  a  reasonably  large  damping  ratio.  Thus,  in  our  attempt  to  create  a  tractable  model,  the  flexible  states 
will  be  removed  from  the  simplified  nominal  model.  In  addition,  the  coupling  of  the  elevator  deflection  5e 
to  the  L  and  D  forces  is  an  undesirable  parasitic  effect.  By  including  this  term,  5e  appears  in  lower  order 
derivatives  of  Vt  and  ft,  requiring  additional  dynamic  extension  to  achieve  vector  relative  degree.  A  very 
similar  situation  involving  parasitic  couplings  in  vertical  take-off  and  landing  (VTOL)  aircraft  is  discussed 
by  Hauser  et  al.  in  Ref.  9.  By  ignoring  the  weak  elevator  couplings,  in  particular  setting  to  zero  the  terms 
Cft,  Cft,  and  Cj,  the  second  order  dynamic  extension  of  Se  becomes  unnecessary.  Figure  2  shows  the  effect 
of  neglecting  these  couplings  on  the  curve  fits  for  L  and  D.  While  still  present,  these  effects  are  clearly 
dominated  by  the  terms  depending  on  a. 


(a)  Drag  Force  with  and  without  Se  coupling 


(b)  Lift  Force  with  and  without  Se  coupling 


Figure  2.  Neglecting  the  coupling  of  Se  to  L  and  D  does  not  significantly  alter  the  accuracy  of  the  fit. 


With  the  removal  of  the  four  flexible  states  and  the  parasitic  elevator  couplings,  the  simplified  nominal 
model  has  only  seven  states:  the  five  original  aerodynamic  variables  and  the  two  integrators  appended  to  the 
<f>  dynamics  from  (5).  The  inputs  for  this  version  of  the  model  are  <hc  and  6e.  Symbolic  calculations  reveal 
that  the  resulting  model  has  full  vector  relative  degree  r  =  7.  The  following  section  discusses  the  control 
design  based  on  this  model. 


IV.  Approximate  Feedback  Linearization  Design 


The  feedback  linearization  design  for  the  simplified  nominal  model  was  completed  using  a  generic  feedback 
linearization  code  for  square  MIMO  systems  that  was  developed  using  the  Matlab  symbolic  toolbox.  Since  the 
system  has  full  relative  degree,  this  design  process  proceeds  using  standard  techniques.10,13  For  simplicity, 
denote  the  state  of  the  simplified  nominal  model  as  x  G  M7.  The  system  can  then  be  rewritten  in  normal 
form  using  a  suitable  change  of  coordinates  to  obtain 


z  =  T(x) ,  T(x) 


zi 


Vt 

LfVt 

L)Vt 

h 

Lfh 

Llh 

L)h 


Z2  =  Z$ 

z3  =  L3fVt  +  (L*cL2fV)$c  +  (LSeL2fV)5e 
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Z4  =  ^5 

Z5  —  Zq 

Zq  =  Z7 

z7  =  Ljh  +  (L^cL)h)^c  +  (. LSeLjh)5e . 


Denoting  the  outer  loop  controls  by  uy  and 
linear 


Uh,  the  following  control  law  renders  the  input/output  map 


Se  _ 
A(x) 


A 


.1  Uy  —  -L^V 
uh  -  Ljh 

L$cL)V  LSeL2fV 
L^c-L^h  L$eLjh 


Specifically,  the  system  under  this  control  can  be  rewritten  as 


k(3) 


Uy 
Uh • 


The  outer  loop  linear  controllers  for  these  two  decoupled  integrator  chains  are  designed  using  LQR  with 
integral  augmentation  and  a  model  reference.  The  LQR  design  weights  are  selected  as 

Qvt  =  diag(10, 1,1,1) 

RVt  =  1 

Qh  =  diag(10, 1, 1, 1, 1) 

Rh  =  1. 


To  obtain  a  reference  model  for  each  linear  system  to  track,  the  reference  signals  for  Vt  and  h  are  filtered 
using  fast  low  pass  filters  with  unity  DC  gain.  The  outputs  of  the  reference  model  are  denoted  Rref  and  hre f. 
By  using  a  filter  of  sufficiently  high  order  for  each  signal,  the  derivatives  of  the  reference  signals  used  below 
are  made  available  to  the  outer  loop  controller.  Specifically,  denoting  with  Ky  and  Kh  the  LQR  design 
gains,  the  outer  loop  controller  is  given  by 


uv{t) 


Uh{t) 


T/(3)  _ 

^ref  ^  ^ 


V  -  Fref 
LfV-Vret 
LjV  -  Kref 

Jo(V(r)  -  Kef (r))dr  . 


1(4) 

*ref 


Kh 


h  href 
Ljh  hre  f 
Ljh  hre  f 
L)h-h^ 
fo(h(T)  ~  Kef(V))dT 


Figure  3  provides  a  conceptual  diagram  of  the  completed  controller. 


Vd 


Figure  3.  The  inner  loop  controller  represents  the  approximate  feedback  linearization  controller  design  based 
on  the  simplified  nominal  model,  while  the  robust  outer  loop  controller  represents  the  outer  loop  LQR. 

The  decoupling  matrix  A(x)  was  verified  to  be  non-singular  over  the  operating  range  of  interest.  Symbolic 
calculations  reveal  that  there  are  four  real  roots  of  the  determinant  of  A.  Two  of  the  roots,  V  =  0  and 
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a  =  0  +  correspond  to  non-physical  situations  and  were  encountered  in  other  similar  models.20  One 
additional  singular  point  corresponds  to  a  an  unattainable  value  of  h.  The  final  solution  yields  somewhat 
plausible  values  for  <a,  but  these  values  are  all  outside  the  anticipated  operating  envelope. 

The  controller  is  guaranteed  by  design  to  perform  well  on  the  simplified  nominal  system.  Simulations 
verify  that  the  controller  linearizes  the  simplified  nominal  model  up  to  the  numerical  accuracy  of  the  algorithm 
used  to  integrate  the  nonlinear  differential  equations,  provided  that  the  variation  of  /?*  is  removed.  (Recall 
from  section  II  that  the  slow  variation  of  /?*  with  Vt  and  h  is  neglected  in  the  control  design.)  Even  with 
this  variation  in  place,  the  input/output  map  is  very  nearly  linear.  Simulation  results  on  the  nominal  model 
demonstrate  good  tracking  performance.  These  results  will  be  omitted  both  for  brevity  and  due  to  their 
similarity  to  the  truth  model  results,  which  are  presented  in  the  next  section. 


V.  Truth  Model  Simulation  Results 


The  truth  model  is  implemented  in  Simulink  using  S-functions.2,8  All  simulations  used  fixed  step  Runge- 
Kutta  integration  with  a  step  size  of  0.05  seconds.  Since  the  truth  model  does  not  support  simulation  of 
high  fuel  to  air  ratios,2  the  attainable  thrust  is  limited  and  all  reference  trajectories  are  passed  through  a 
rate  limiter  followed  by  a  low  pass  filter  of  the  form 


H(s) 


.04 

s2  +  As  +  .04 


The  selected  reference  trajectory  begins  at  q  =  2000^  and  h  =  85000  ft.  The  aircraft  climbs  at  a  steady 
50^  while  maintaining  a  constant  dynamic  pressure.  Once  the  HSV  reaches  mach  10,  the  mach  number 
is  held  constant  and  the  climb  rate  increases  to  139^  until  levelling  off  at  115,000  ft.  This  trajectory  was 
determined  to  be  a  plausible  operating  trajectory  for  the  vehicle,  and  the  ramp  increase  in  altitude  provides  a 
moderately  aggressive  tracking  challenge  to  the  controller.  It  should  be  noted  that  the  available  thrust  limits 
the  climb  rate  for  this  reference  trajectory.  If  a  faster  climb  rate  is  selected  for  the  first  leg,  the  required 
acceleration  cannot  be  achieved  without  violating  the  actuator  limits  on  4>.2 

Figures  4  and  5  show  the  simulation  results  for  the  test  trajectory  on  the  truth  model  under  the  proposed 
controller.  The  tracking  performance  is  excellent  throughout  the  entire  maneuver,  and  both  the  flexible 
states  and  pitch  dynamics  remain  stable  and  well-behaved.  Simulations  were  also  conducted  with  mild  plant 
parameter  variations  to  quantitatively  assess  the  robustness  of  the  controller.  This  brief  case  study  revealed 
that  the  controller  is  at  least  moderately  robust.  Figure  6  gives  an  example  result  from  this  study  in  which 
the  vehicle  mass  was  increased  by  roughly  10%,  Iyy  was  decreased  by  25%,  and  the  vehicle  length  was  reduced 
by  15%.  The  tracking  performance  is  still  very  good,  the  most  noticeable  difference  being  an  overshoot  in 
Vt  when  the  trajectory  reaches  mach  10  and  levels  off. 


VI.  Control  Design  for  an  Alternative  Output  Choice 

While  the  output  choice  Vt  and  h  is  quite  logical  for  an  auto-pilot,  a  controller  designed  for  Vt  and  the 
flight  path  angle  7  would  be  more  logical  for  use  by  a  pilot.  With  this  output  selection,  the  simplified 
nominal  model  has  vector  relative  degree  r  =  6,  the  one  dimensional  internal  dynamics  becoming  that  of  the 
vehicle  altitude  h.  However,  on  the  output-zeroing  manifold  7  =  0,  and  examination  of  equation  2  reveals 
that  the  zero  dynamics  are  given  by  the  trivial  equation  h  =  0.  Thus  the  altitude  will  remain  stable  even 
though  it  is  rendered  unobservable  by  the  linearizing  feedback.  Figure  7  and  8  show  simulation  results  for 
the  controller  using  Vt  and  7  as  the  regulated  output.  As  expected,  the  trajectory  of  the  altitude  h  remains 
bounded.  It  should  be  noted  that  the  test  trajectory  used  for  this  controller  is  a  simple  step  increase  in  Vt 
and  a  unitary  pulse  for  7.  Also,  the  outer  loop  controller  was  implemented  using  a  simple  pole  placement 
design.  A  more  sophisticated  outer  loop  design  similar  to  the  one  employed  for  the  Vt  and  h  output  choice 
could  yield  superior  performance  for  this  choice  of  output. 
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Altitude,  ft 


(a)  Altitude 


(b)  Velocity  Result 


(C)  Se  (d)  4> 

Figure  4.  Simulation  Result  for  Baseline  Trajectory  on  the  Truth  Model 
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(a)  Angular  States 


(b)  Flight  Path  Angle 


-  drij/dt  A  x1 ,  deg/sec 

-  dr|  /dt  A  t„,  deg/sec 


So- 
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f 

0  50  100  150  200  250  300  350  400  450 

Time,  sec 


(c)  Flexible  States  (d)  Flexible  State  Rates 

Figure  5.  Simulation  Result  for  Baseline  Trajectory  on  the  Truth  Model 
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Altitude,  ft 


(a)  Altitude 


(b)  Velocity  Result 


(C)  Se  (d)  $ 

Figure  6.  Simulation  Result  for  Baseline  Trajectory  on  the  Truth  Model  with  Parameter  Variation. 
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(a)  Flight  Path  Angle  (b)  Velocity  Result 


(c)  Se  (d)  <F 

Figure  7.  Simulation  Result  for  7  and  V  Tracking 
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(a)  Angular  States 
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(c)  Flexible  States  (d)  Flexible  State  Rates 

Figure  8.  Simulation  Result  for  7  and  V  Tracking 
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VII.  Conclusions 


In  this  paper  we  have  demonstrated  the  suitability  of  control  system  design  based  on  approximate  feedback 
linearization  for  an  air-breathing  hypersonic  vehicle.  The  key  steps  in  the  design  are  the  construction  of  an 
analytically  tractable  nominal  model  and  the  removal  of  certain  flexible  effects  and  weak  couplings  to  achieve 
full  relative  degree  for  a  simplified  model.  The  resulting  controller  achieves  comparable,  if  not  superior, 
performance  to  earlier  linear  control  designs8  over  a  wider  operating  range.  In  addition,  since  this  approach 
is  more  closely  tied  to  the  nonlinear  model,  an  extensive  tuning  algorithm  was  not  necessary  to  achieve 
excellent  tracking  performance.  Of  course,  the  penalty  for  this  more  straightforward  outer  loop  design  lies 
in  the  dramatically  increased  complexity  of  both  the  system  identification  process  and  the  computational 
complexity  of  the  controller  itself. 

The  simplified  nominal  model  was  constructed  in  the  spirit  of  approximate  feedback  linearization  as 
proposed  by  Hauser  et  al.9  Future  work  will  include  an  analytical  proof  of  stability  for  the  nominal  system 
under  the  designed  control  law,  along  with  an  attempt  to  generalize  this  result  to  a  class  of  similar  systems. 
In  addition,  more  elaborate  vehicle  models  with  additional  flexibility  effects,  thermal  effects,  and  engine 
dynamics  will  be  considered.  Additional  control  effectors,  such  as  the  canard  used  in  a  recent  paper  by 
Bolender,7  may  also  be  employed  to  improve  control  performance  for  these  more  sophisticated  models. 
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